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We report on a systematic investigation of the critical fluctuations in the complex conductivity 
of epitaxially-grown Bi2Sr2CaCu208+i films for T <T C using a two-coil inductive technique at zero 
applied field. We observe the static 3D-XY critical exponent in the superfluid density near T c . Linear 
scaling analysis close to the critical temperature yields a dynamic critical exponent of z ~ 2.0 for 
small drive currents, but non-linear effects are seen to be important. At T c , a non- linear scaling 
analysis also yields a 3D dynamic exponent of z = 2.0 ± 0.1. 
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I. INTRODUCTION 

Direct observation of the scaling properties of the 
superconducting transition has at last become feasi- 
ble due to the discovery of the cuprate superconduc- 
tors. Microwave measurements in ultrapure single crystal 
YBa2Cu307_a; (YBCO) have provided strong evidence 
that for accessible temperature ranges, the effective static 
universality class is the three dimensional XY model (3D- 
XY)A>, where the penetration depth displayed scaling over 
three decades in reduced temperature. Measurements 
on thin films have yielded varying results: microwave 
measurements* report v = 1.2, contactless ac conductiv- 
ity measurements* yield v = 1.7 at frequencies up to 2 
GHz and dc conductivity in finite magnetic fields* give 
v = 0.9 — 1.0. While none of these agree with the 3D-XY 
scaling behavior, it is widely believed that the results on 
pure single crystals are indicative of the intrinsic critical 
fluctuations of the superconducting transition. 

Transport measurements have probed the dynamics of 
the critical fluctuations, which are characterized by the 
value of the dynamic critical exponent z describing how 
the relaxation time r scales with the correlation length 
£: t ~ £ z . However for YBCO, experiments have not 
yielded a consistent picture. For example, longitudi- 
nal dc-resistivity measurements* yield z — 1.5±0.1 in 
magnetic fields up to 6T, while microwave conductiv- 
ity measurements* are consistent with z — 2.3 — 3.0. 
Still larger values for the dynamic exponent ranging from 
z = 5.6 to z = 8.3, reminiscent of glassy dynamics, have 
also been seen in resistivity measurements** in thin films. 
These large values of z strongly suggest that disorder 
plays a dominant role in these thin films, and this may 
explain why their observed static critical behavior differs 
from that of the best single crystals. 

Critical fluctuation results on Bi 2 Sr2CaCu208+,5 
(BSCCO) are even more dissimilar, partially due to 



highly anisotropic fluctuations. The complex con- 
ductivity of underdoped BSCCO films measured at 
terahertz frequencies was found to agree with a 
Kosterlitz-Thouless-Berezinskii (KTB) model with diffu- 
sive dynamics!** 1 In previous work on BSCCO crystals 
a non-linear transport study found 2D critical behavior 
with an anomalous dynamic exponent of z ss 5.6 and no 
crossover to 3D behaviori* The 3D universal phase an- 
gle of the complex conductivity deduced from the mag- 
netic susceptibility of BSCCO epitaxial films and crys- 
tals has even yielded larger values of zi Only dc fluc- 
tuation conductivity measurements in BSCCO** find re- 
laxational dynamics, z « 2.0, assuming a 3D-XY static 
exponent, v w 2/3, in agreement with recent theoretical 
work that predicts that z s» 2.0 in 3D *i* These results in- 
dicate the lack of consensus on the true nature of the dy- 
namic fluctuations representative of the superconducting 
transition, which is the primary motivation of the present 
investigation. While a large value of the dynamic expo- 
nent (z > 2) can be understood as disorder effects, only 
the presence of a previously undetected collective mode 
coupled to the superfluid density would lead to z ss 1.5 
in 3D, just as occurs in superfluid He 4 >i* 

In this paper we provide results on the critical fluctua- 
tions in high quality epitaxially-grown BSCCO films. We 
have measured the ab-plane complex conductivity with 
a two-coil inductive technique at zero applied magnetic 
field (H = 0), described in Section II. Our measurements, 
presented in section III, indicate that the critical fluctu- 
ations in BSCCO films are consistent with the 3D-XY 
critical fluctuation model, rather than a KTB transition 
within the layers. This contrasts with earlier studies us- 
ing this technique on YBCO, which have reported KTB 
critical!* and mean-field** fluctuations. In section IV, 
we estimate the temperature intervals in which possible 
crossovers to 2D fluctuations can occur, due either to de- 
coupling of the bilayers or to the thin film behaving as a 
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2D system with a large c-axis correlation length. We find 
that our measurements should indeed be well-described 
by anisotropic 3D-XY critical fluctuations. In section V, 
we provide a crossover analysis of the approach to the 
critical point, and examine the approach to a universal 
phase angle. In particular, we predict the qualitative 
forms of the variations with temperature and frequency. 
These are in agreement with our measurements, but con- 
trast earlier results which report very different values of 
the static and dynamic critical exponents than those re- 
ported here. Analysis of the phase angle of the film re- 
sponse using linear response theory in the critical region, 
which yields a dynamic exponent of z = 2.0, is presented 
in section VI. However, the breakdown of linear response 
theory becomes important near the transition tempera- 
ture T c , so that the analysis is complicated by the need to 
extract the superfluid density from the measured mutual 
inductance by an approximate inversion technique. Our 
analysis shows that the determination of z purely from 
linear response theory is contaminated by non-linear ef- 
fects, and at probe currents that are too large the effec- 
tive value of z inferred can be as low as 1.5. To circum- 
vent this problem, we derive in section VII a non-linear 
scaling law obeyed by the raw measured mutual induc- 
tance at T c (i.e. without the need for a data inversion) 
and extract the dynamic exponent directly from the mu- 
tual inductance. In agreement with the linear response 
analysis, we again obtain z ss 2.0. Section VIII summa- 
rizes our conclusions. 



II. SAMPLES AND EXPERIMENTAL SETUP 

We have measured several high quality oxygen-doped 
BSCCO films, grown by molecular-beam-epitaxy on 
SrTiC>3 substrates and analyzed in- situ by RHEED 16 . 
The BSCCO films are nominally oriented along the c- 
axis, with unit cell thickness of d = 15.4 A. AFM images 
of similar BSCCO films reveal an rms roughness of 5 
Angstroms primarily due to single unit cell step height 
variation over the surface. The rms roughness of the 
substrates are approximately 1.5 Angstroms. 

The inductively-measured T c from several films of dif- 
ferent dopings was fit to the empirical formula T c — 
T c ,max{^ — 82.6(p — 0.16)) 2 , to determine the hole doping 
per Cu, p^ The two branches of the formula are distin- 
guished by the relative low temperature superfluid den- 
sity, which is known to increase monotonically through 
optimal doping (p=0.16). Films A, B, and C have 
ii = 21,40 and 60 unit cells respectively and are grown 
on 10 and 14 mm square substrates. Film A is overdoped 
at p ~ 0.19, while film B and C are slightly overdoped 
and underdoped, respectively. We estimate that a dop- 
ing gradient in the samples introduces a experimental 
temperature broadening of 0.15 K near T c . 

To obtain the ab-plane complex conductivity, we em- 
ploy a two-coil inductive technique. The mutual induc- 
tance of two coaxial coils is measured with a film placed 



in between and normal to the axis of the coils. An ac 
current is applied through the drive coil; a second coil, 
attached to a large input impedance lock- in amplifier, 
measures the fields produced by the drive coil and the 
screening currents in the film. The coils have 135 turns 
and an average radius of 1.5 mm. The drive coil is placed 
0.35 mm above the film so that at a current of I = 50/iA 
rms, the applied ac magnetic field is < 0.1 Gauss nor- 
mal to the film. Measurements are performed in a He 4 
continuous- flow dewar, which allows the temperature of 
the film to be controlled without significant heating from 
the drive coil. The sample temperature is monitored by 
a Si diode attached to the back side of the substrate. 

The complex conductivity a = <j\ — i<72 of each film 
is calculated from the mutual inductance using the ex- 
act analytical expression for an infinite diameter film 17 
and a numerical inversion technique similar to that of 
Fiory et ali&. The method of Turneaure et nlJ^ is em- 
ployed to correct for film diameter effects prior to calcu- 
lating the in-plane complex conductivity. The measured 
penetration depth A a f, = (po UJ &2)~ 1 ^ 2 , is independent of 
the measurement frequency from / = 10 — 100 kHz be- 
low the transition temperature, where the dissipation is 
small (a i <C 02). To analyze the critical fluctuations, we 
study the complex superfluid density per CuO bilayer, 
p = ioujd§\l (47r 2 fcs), where <I>o is the superconducting 
flux quantum^ The real part of this quantity, i£e[p] has 
units of temperature and sets the energy scale for critical 
fluctuations. In Fig. 1(a), Re[p] obtained at / = 80 kHz 
and I = AOpA is plotted for films A through C. At 5 K, 
films A, B, and C, have penetration depths of 235, 265, 
and 285 nm, respectively. Note that film A, which is over- 
doped, has a smaller low temperature penetration depth 
than film B and C. All films exhibit a T 2 component to 
the superfluid density at low temperatures, however the 
zero-temperature intercept of the films near optimal dop- 
ing is approximately zero, indicating the crystalline films 
have very few impurities. 



III. STATIC CRITICAL BEHAVIOR 

In Figs. 1(a) and 1(b) a power of p = Re[p] and 2T/7T is 
shown to locate where we would expect the KTB transi- 
tion for a noninteracting bilayer, p = 2T/ir*Q We observe 
no such drop in the superfluid density at this tempera- 
ture in our films. The films instead show that coupling 
among the bilayers induces an anisotropic 3D transition, 
precluding a KTB transition in the bilayers. At the tran- 
sition temperature of an isolated bilayer T — pir/2, how- 
ever, the measured superfluid density has been renormal- 
ized from the mean field value by quasi-2D fluctuations. 

The mean field behavior, expected at low tempera- 
tures, crosses over to the 3D-XY behavior as the crit- 
ical temperature is approached from below. In layered 
systems, the 3D-XY behavior is exhibited once the c- 
axis correlation length exceeds the interlayer spacing, d; 
i.e. when p(T) = c s T c , where c s m 0.5 and p — Re[p\. 1 
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v. At higher temperatures where we observe non-linear 
behavior as well as frequency dependence, the data are 
still consistent with the 3D-XY fluctuations, and we can 
extract z. In section VI we analyze the phase angle in 
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FIG. 1: Re[p] (panel a) and Re[pf /2 (panel b), as a function 
of temperature for films A through C. p = (c s 85K) is the pre- 
dicted crossover to 3D-XY critical behavior and p = (2/n)T 
is the KTB critical superfluid density for an isolated bilayer. 



In Fig. 1(b), the data from Fig. 1(a) are plotted as 
i?e[p] 3 / 2 near the critical regime. Also shown in Figs. 
1(a) and 1(b) is the corresponding power of c s T c to lo- 
cate the onset of 3D-XY behavior. The static 3D-XY 
model gives p ~ £ _1 ~ |T - T c \ v , where v « 2/3. We 
observe that Re[p] 3 ^ 2 varies linearly in temperature for 
p(T) < c s T c indicating a 3D-XY, rather than mean-field 
(p ~ \T — T c \), behavior. At our measurement current 
(40/iA), we observed an absence of frequency dependence 
below the dissipation peak over the temperature range of 
the fit. In film C, the linear regime extends as far as 7 K 
below 82 K. The linear fits extrapolate to T c = 79.1 K , 
84.9 K and 84.7 K for films A, B and C respectively. 

Although our static scaling analysis is carried out at 
temperatures and at a drive current for which there is 
no discernible frequency dependence, at higher tempera- 
tures we observe variations with both drive current and 
frequency. This is apparent in Fig. 2, for which larger 
currents are used on film B with 40 layers. Film C, with 
60 layers, exhibited a weaker dependence over the same 
frequencies and currents. Since the superfluid density de- 
pends on frequency and current at higher temperatures, 
the superfluid density scales as p <~ £ 2 ~ D f(oj£ z , E£ 1+z ), 
where E is the amplitude of the electric field. Only in the 
critical regime where the data are independent of both 
frequency and electric field can one expect to see the lin- 
ear response scaling behavior and extract the exponent 




FIG. 2: \p\ in film B taken at a series of frequencies (panel a) 
and drive currents (panel b). 



IV. DIMENSIONALITY OF THE STATIC 
FLUCTUATIONS 

When the superfluid density is a fraction of the 3D 
crossover (p(T) < (l/2)c s T c ), we observe a drop in the 
superfluid density from the static scaling value. Although 
this could be interpreted as a crossover to 2D behavior, 
we do not believe this to be the case in our films since the 
drop is observed to be current as well as frequency depen- 
dent (see Fig. 2). Notice that the superfluid density does 
drop more rapidly than expected from the 3D-XY model 
as the frequency is decreased for a given drive current, 
but does so more slowly as the current is decreased at 
fixed frequency While the former might lead one to con- 
clude that the anisotropic 3D-XY behavior is no longer 
valid in the static limit, the latter would imply the op- 
posite. More accurate measurements are needed at low 



4 



frequency and currents to establish unambiguously the 
true nature of the static limit. Nevertheless the trend 
with decreasing current, coupled with the phase angle 
dependence discussed in the next section, is consistent 
with 3D-XY behavior. Another putative crossover to 2D 
behavior is expected when £ c becomes as large as the film 
thickness, h. The 3D-XY scaling implies that this occurs 
at a temperature, T*, which is 0.05 K, 0.01 K and 0.02 
K below the transition temperature for films A, B and 
C respectively. This implies that the region of pure 2D 
fluctuations is too small to be resolved in our data and 
therefore we expect 3D fluctuations to dominate. 
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FIG. 3: Normalized phase angle, 2<j> p /n, of the superfluid 
density extracted from linear response theory. Panel a: 2<f> p /n 
in film B taken at / = 100/iA (filled symbols) at a series of 
frequencies and I — 50/j, A for / = 50kHz (unfilled triangles). 
Panel b: 24> p /n in film C for a given table of frequencies and 
drive currents. 



V. CROSSOVER PHENOMENA IN THE 
LINEAR PHASE ANGLE NEAR T c 

Since the experiments are performed at finite fre- 
quency, it is important to consider whether or not the 
correct asymptotic behavior is being probed and what 



conclusions can be drawn from the observed behavior. 
We begin with the superfluid density, within linear re- 
sponse, written as 

P = (t-»Y- D h{ut-» Z ) (i) 

where t — \T — T c \ /T c and /i is a scaling function. In 
the limit of t — > (T — * T c ), the i-dependence must 
cancel out of the expression for p, which implies that the 
asymptotic behavior of the scaling function fi(x) ~ x e 
as x — > oo where 9 = (2 — D)/z. For t 0, we can write 

p={uj)P- 2 V'f2{tu- 1 ' v *) (2) 

where the scaling function = const as y — > 0. We 

now ask what behavior we should observe as ui — ► for 
i ^ 0. For oj — > at t — 0, we should see the criti- 
cal point behavior. For t/uj 1 ^ 2 C 1 we are probing the 
y ~ regime of the scaling function f2(y) and expect 
to see asymptotic scaling behavior. If t/ui 1 / vz » 1 wc 
are probing the y — -» oo regime of the argument of /2 
and do not expect to see the scaling behavior. Thus at 
a finite to, for uj 3> t^ 3 we would observe critical be- 

4/3 

havior while for u> <C t we would not. We have used 
the critical exponents v — 2/3 and z = 2 for specificity, 
but the argument does not rely on these precise values. 
Note that, as is generally the case with crossover, the 
behavior of the superfluid density away from the critical 
point at finite frequencies is not universal except in the 
large frequency limit - a result that is sometimes found to 
be surprising, but which follows from a generic crossover 
analysis of relevant variables. 22 This predicted low fre- 
quency behavior is consistent with our data where we 
observe the critical phase angle in Fig. 3(b) and yet we 
do not observe scaling below T c . More importantly, this 
is also consistent with the non-linear scaling observed at 
T c in Section VII. Away from T c an asymptotic scaling 
regime is demonstrated in the microwave measurements 
of Booth et al& 



VI. DYNAMIC CRITICAL BEHAVIOR AND 
THE BREAKDOWN OF LINEAR RESPONSE 
THEORY 

Next we study the phase of p, <f> p , close to the tran- 
sition temperature. The frequency dependence of 24> p /ir 
for film B measured at / =100 pA is shown in Fig. 3(a). 
Within linear response theory, the phase angle is inde- 
pendent of frequency at the critical temperature and is 
given by <p p (T c ) = n/(2z) in three dimensions, where 
z is the dynamic critical exponent^ In Fig. 3(a) we 
notice that the curves for different frequencies seem to 
approach a phase angle consistent with a dynamic expo- 
nent of z w 1.5 and T c = 84.9 K. When repeated at a 
smaller value of the current, I = 50pA, 4> P (T C ) is smaller, 
as shown. Since the constant phase angle is a result from 
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linear response theory, the estimate of the dynamic ex- 
ponent from the data is expected to be more accurate at 
smaller currents. 

This is indeed seen in the current dependence of the 
phase angle at different frequencies for film C. In Fig. 
3(b), 2<p p /w is shown for a set of frequencies and cur- 
rents. As the current is lowered, the frequency inde- 
pendent phase angle shifts to lower values. At 20/iA, 
lOpA and 5/zA, this phase angle is 0.577r/2, 0.547r/2 and 
0.5i7r/2 respectively. This suggests that z ~ 2.0 in film C 
and the result of z » 1.5 from film B is an artifact of the 
large current used in the measurement. With this anal- 
ysis, the thickest film (C) exhibits the best estimate of 
the critical phase angle, the film with intermediate thick- 
ness (B) provides evidence for similar behavior, but the 
thinnest film did not show a critical phase angle. Due to 
the crossover phenomena discussed in the previous sec- 
tion, the behavior of the phase angle away from T c is not 
universal at low frequencies, but we can still infer the 
dynamic critical exponent from the critical phase angle 
at T c which is indeed universal. 

Both the phase angle and the superfluid density are 
found to be dependent on the drive current near T c . In 
the limit of smaller currents, the frequency independent 
critical phase angle does indeed yield a dynamic exponent 
of z Kt 2.0. This is a result consistent with the presence of 
3D fluctuations. This is also consistent with the behavior 
of the superfluid density with decreasing current at fixed 
frequency, where it is seen to approach the asymptotic 
static scaling law. Furthermore, the smooth evolution 
of the phase angle suggests that the crossover to the 2D 
regime near T* is not resolved in our experiment. Nev- 
ertheless, considering that nonlinear scaling and dimen- 
sional crossover is expected close to T c , we now proceed 
to analyze the nonlinear response to get a more accurate 
determination of the dynamic exponent. 

VII. NONLINEAR SCALING AT THE 
SUPERCONDUCTING TRANSITION 

The field equation governing the vector potential, A, 

is 

V 2 A//i = -J d - hS(z)aE (3) 

where Jd is the current density in the drive coil, a is the 
conductivity of the superconducting film, and E is the 
electric field. For the geometry of the setup, 

Jd = IdS{r - r c )S(z - z c )$ (4) 

in cylindrical coordinates, where Id is the drive current, 
r c is the radius of the drive coil and z c is the distance from 
the film to the drive coil^i For an anisotropic system, in 
the critical regime, the superfluid density, p, scales as, 

p ~ iua ~ e^^f^Qb, E£t') (5) 



where and £ c are the correlation lengths parallel and 
perpendicular to the CuO bilayers, and / is a scaling 
function. We can rewrite the field equation as 

V 2 A = - p I d S(r - r c )S(z - z c ) (6) 

where we have used the scaling form for the complex con- 
ductivity. Since the scaling function is written in terms 
of dimcnsionless variables we can non dimensionalize the 
equation by recasting it in terms of A£ a b, u)£* b , r/h and 
Id£ab- The solution then takes the form, 

= G{iu£ b ,iue a i z h,h%^ D Z- 1 A.-\r/h) (7) 

where G is a function that can be determined by solving 
the full non-linear field equation and A is the thermal 
length given by $§/(47r/iofcsT). So far all we have done 
is to use the scaling form for the superfluid density and 
looked for a generic solution of a differential equation 
by identifying all the non-dimensional variables. For a 
given geometry of the experiment the right hand side of 
Eq. [7| is a function of three variables. One can further 
simplify the expression if we realize that there exists a 
regime \T — T c \ < T* where we are governed by the limit 
where £ c = h. This allows us to eliminate one of the three 
terms as follows. Close to T c , using the 3D-XY behavior 
of the correlation length and the definition of T*, we can 
rewrite h 2 ^ D /^c as, 

h 2 Cb°/^ = h\l + (T — T*)/(T* - T c )\ v . (8) 

For a temperature range \T - T*\ <C \T* - T c \, h 2 /^ c « 
h. Thus, for temperatures within 0.05AT, Q.02K, and 
0.01K of T* for films a, b, and c respectively, this ap- 
proximation is valid, and the field in Eq. 0is a function 
of only two variables for a given geometry. Thus one 
can now perform the standard data collapse at T c of the 
measured mutual inductance to obtain the dynamic ex- 
ponent. In this regime we can eliminate t; a b from Eqn0 
so that 

A = uj 1 / z G(I d oj- 1/z ,hA~ 1 ,r/h). (9) 

The mutual inductance, M, is related to the vector po- 
tential at the pick up coil and scales as 

M ~ Alj- 1/z ~ G{I d iuj- 1/z ). (10) 

Within the resolution of our data, any data collapse 
observed is a 3D phenomena starting to crossover to the 
asymptotic 2D scaling. The reason for this is that the 
solution (Eq. UJ depends on a number of variables, but 
sufficiently close to T c there exists a regime wherein the 
mutual inductance is only a function of one scaling vari- 
able. If we do not observe any data collapse we would 
conclude that we are never in the regime where the ap- 
proximation h 2 /£ c Ri h is valid. It is indeed true that 
the elimination of the dependence of £ a j, in Eq. {7\ is pos- 
sible only at T c and need not hold in the entire regime 
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where the crossover occurs. It is tempting to interpret 
the experimental results in terms of a 2D system. Our 
data on the other hand are unable to resolve the crossover 
regime, and yet do exhibit data collapse (see Fig 4). This 
suggests that the temperature at which we observe data 
collapse lies in the crossover regime, but is not necessarily 
the true T c . From our analysis we conclude that the T c 
is within 0.02 K of the temperature at which we observe 
data collapse for film B and C. We now look for data col- 
lapse of our data measured by varying frequency, current 
and temperature. Notice that the value of the exponent 
v is required only to establish the regime of validity of 
the approximation and is not necessary to extract the 
dynamic exponent from data collapse at T c . 
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FIG. 4: Panel (a): \M\ measured at I = 100^A and / = 50 
kHz versus the scaling variable Icu^ 1 ' 2 . Panel (b): Scaling 
error at different temperatures versus z. 

For the scaling analysis at T c , the mutual inductance 
at fixed frequency, M/, was compared to the mutual in- 
ductance at fixed current, Mj. For convenience, \M(T)\ 
is normalized to unity above the T c . Values of Mf (Mi), 
were measured at 12.5 /iA (10 kHz) intervals for a set 
of fixed temperatures. Then values of \Mj\ were com- 
pared with values of |M/| taken at the same temperature 



by fitting the raw data of \Mf\ to a polynomial and se- 
lecting points with the values of Ioj^ 1 / 2 equivalent to 
those of \Mi\. In Fig. 4(a) \M f \ and |M 7 | for film B are 
plotted versus Ioj-V 2 for a series of temperatures. The 
best data collapse of the curves appear at z — 2.0 and 
T = 85.22 K. The error in the scaling, Am, is shown 
in Fig. 4(b) for different temperatures and values of z. 
This is obtained from A M = J2(\ M f \ ~ \ M i\) 2 , where 
the sum runs over the 10 values of Iuj^ 1 ^ measured for 
\Mj\. The lowest error value indicates z — 2.0 ± 0.1 and 
T = 85.2 ± 0.1 K. A similar analysis for the film C yields 
z = 1.8±0.2 and T = 85.9±0.2 K, which is also consistent 
with z — 2.0 ±0.1. The transition temperature of film B 
is T c = 85.2 ± 0.1K, while for film C it is T c = 85.9 ± 0.2 
K. These values of T c are close to those obtained with 
the phase angle measurements. 

In Fig. 4(b), the evolution of the error for a fixed trial 
z also suggests that the dynamic exponent observed is 
indeed due to 31? fluctuations. Notice that as the tem- 
perature is increased for z = 2.0, the error decreases, 
approaching the best fit (i.e. data collapse), and then 
increases again. Within the resolution of the data, we 
never observe any change to the 2D regime. In other 
words, the approximation that allowed us to look for the 
data collapse is indeed supported by the behavior of the 
mutual inductance. 



VIII. CONCLUSIONS 

The superfluid density in BSCCO films with different 
thicknesses has been measured and shown to be consis- 
tent with the predictions of the 3D-XY model. The dy- 
namic scaling exponent has been obtained by performing 
both a linear and non-linear scaling analysis. While the 
linear scaling analysis did yield a 3D dynamic exponent 
of z w 2.0, non-linear effects needed to be studied to 
get an accurate determination. This is evident from the 
current dependence of the critical phase angle of the su- 
perfluid density. A 3D non-linear scaling analysis also 
yields a dynamic exponent of z = 2.0 ± 0.1. 
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